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ABSTRACT 

The recent emergence of a new class of accretion-powered, transient, millisecond X-ray pulsars presents 
some difficulties for the conventional picture of accretion onto rapidly rotating magnetized neutron 
stars and their spin behavior during outbursts. In particular, it is unclear from the standard paradigm 
how these systems manage to accrete over such a wide range in M (i.e., > a factor of 50), and why 
the neutron stars exhibit a high rate of spindown in at least a number of cases. Following up on prior 
suggestions, we propose that 'fast' X-ray pulsars can continue to accrete, and that their accretion 
, disks terminate at approximately the corotation radius. We demonstrate the existence of such disk 

■ solutions by modifying the Shakura-Sunyaev equations with a simple magnetic torque prescription. 

The solutions are completely analytic, and have the same dependence on M and a (the viscosity 
parameter) as the original Shakura-Sunyaev solutions; but, the radial profiles can be considerably 
modified, depending on the degree of fastness. We apply these results to compute the torques expected 
during the outbursts of the transient millisecond pulsars, and find that we can explain the large 
spindown rates that are observed for quite plausible surface magnetic fields of ~ 10 9 G. 
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1. INTRODUCTION 



The recent discovery of five accretion-powered millisecond X-ray pulsars (SAX J1808-3658; XTE J0929-3I4; XTE 
q JI751-305; XTE 1807-294; XTE J1814-338; Wijnands & van der Klis 1998; Chakrabarty & Morgan 1998; Galloway et 
■ al. 2002; Markwardt et al. 2002; Markwardt, Smith, & Swank 2003; Markwardt, Juda, & Swank 2003; Markwardt, & 
Swank 2003; Markwardt, Strohmayer, & Swank 2003), all of which are X-ray transients, raises a number of interesting 
C$ • questions about the fundamentals of accretion onto rapidly rotating neutron stars. These 5 msec pulsars have typical 
luminosities that range from about ~ 5 x 10 36 ergs s _1 near the start of the X-ray outbursts, down to as low as 
. £h ! ~ 5 x 10 34 ergs s _1 (for assumed distances of ~5 kpc) as they fade in intensity below the detection limit (e.g., of the 
RXTE satellite). The X-ray pulsations are typically detected as long as any DC intensity is detected from the source. 
' The pulse profiles are basically sinusoidal and do not change dramatically in character as the source intensity drops. 
C3 In particular, there is no evidence of anything like a "propeller effect" setting in as the intensity gradually falls to quite 
low values. 

These observational facts raise a number of puzzling issues (see also, e.g., Burderi & King 1998; Psaltis & Chakrabarty 
1999). First, how are these rapidly rotating, magnetized neutron stars able to accrete matter over such a wide range 
(e.g., ~ 50) in mass accretion rate, M, without encountering a "centrifugal barrier" effect (e.g., Illarionov & Sunyaev 
1975, hereafter, IS; see also Davidson & Ostriker 1973; Fabian 1975; Lipunov & Shakura 1976)? Second, why do at 
least two of the msec X-ray pulsars exhibit a large spindown in pulse period, i.e., P > 0, as the luminosity declines 
throughout the outburst (Galloway et al. 2002; Morgan, Galloway, & Chakrabarty 2003)? Third, why is it that — thus 
far — only relatively faint X-ray transients have been found to exhibit coherent periodic msec pulsations, while many 
other "steady" low-mass X-ray binaries (LMXBs) do not? 

In the conventional picture of accretion onto a magnetized neutron star via an accretion disk, the magnetosphcric 
radius is taken to be 



where M is the mass of the neutron star, M is the steady-state mass accretion rate, and \i is the magnetic dipolc 
moment of the neutron star (e.g., Lamb, Pethick, & Pines 1973 "LPP"; Rappaport & Joss 1976; Ghosh & Lamb 1979 
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"GL"; Wang 1987). The quantity /Z26.5 is the dipole moment expressed in units of 10 26 ' 5 G cm 3 , which corresponds 
to a surface magnetic field at the poles of <~ 3 x 10 s G. Here and throughout this work we take eq. (1) to be a formal 
definition of r m . Inside the magnetosphcric radius the dynamics of the accreting matter are supposedly dominated by 
the magnetic field of the neutron star. The corotation radius is defined as the radial distance at which the Keplerian 
angular frequency is equal to the spin frequency of the neutron star: 

/GM\ 1/3 ol , / P \ 2/3 / M \ 1/3 



According to conventional wisdom, the magnetosphcric radius, r m , is supposed to lie within the corotation radius, 
r c (LPP; GL), in order for accretion to take place. Otherwise, it is thought that the centrifugal barrier, experienced by 
particles forced to corotate faster than the Keplerian velocity, would expel matter from the system (IS). Such pulsars 
are termed "fast" , and arise, for a given system, when the accretion rate drops below the value implied by setting 
r c = r m . At the other extreme, if the accretion rate becomes too large, then the value of r m could drop below a value 
of <~ 10 km, the approximate radius of the accreting neutron star. If the accretion disk persists down to the surface 
of the neutron star, detectable pulsations might not arise. Thus, for typical system parameter values associated with 
the msec X-ray pulsars, r m might have to lie within a rather restricted range, e.g., 10 — 35 km. A range of 3 — 4 in r m 
corresponds to a range in M of about a factor of 100 (see eq. [1]). However, the observed range in X-ray luminosity, 
and presumably also in M, is of order a factor of <~ 50 for at least three of the msec pulsars. Thus, in the standard 
accretion paradigm the magnetic field and values of M would have to be in just the correct range to allow accretion to 
continue throughout the transient outburst (see also Burderi & King 1998; Psaltis & Chakrabarty 1999). Of course, 
this is still plausible as the result of observational selection effects, i.e., the sources we observe to pulse have just the 
right parameters; however, the range of allowable parameter space is growing uncomfortably small. Moreover, the 
neutron star should be observed to be spinning up during most of the outburst, while the opposite was true on at least 
one occasion each, for two of the sources (SAX J1808-3658 and XTE J0929-314). 

The conventional "propeller" picture, in which all mass is expelled as soon as the nominal magnctosphere radius 
lies outside the corotation radius (IS), is known to be incorrect from a theoretical point of view. Spruit and Taam 
(1993, hereafter ST) have pointed out that the velocity excess of the stellar rotation over the Keplerian velocity at r m 
is energetically insufficient to gravitationally unbind all accreting matter, unless r m /r c is actually rather large. Hence 
there must be a regime with r m > r c where accretion in the standard way is not possible, but at the same time the 
bulk of the accreting mass cannot leave the system either. ST showed, with a time-dependent calculation, that what 
happens instead is that accretion continues at some level, while outside r c mass builds up to high surface densities 
(more accurately: surface density normalized to accretion rate). The high densities in the inner disk cause the net 
angular momentum flux to be directed outward, while mass accretion continues inward. 

In this paper we argue that the accretion disk structure around a fast pulsar will adjust itself so that the inner edge 
of the disk will remain fixed near ~ r c , and accretion will continue. Matter will then not be ejected from the system 
in substantial quantity by the propeller effect (IS; if the magnetic moment of the neutron star is not perfectly aligned 
with its rotation axis), in conjunction with a "centrifugal barrier". We compute approximate disk solutions that allow 
for penetration of the disk to r c . The solutions for the pressure, density, temperature, and disk thickness are in the 
form of analytic functions of radial distance from the neutron star and also depend on M, the viscosity parameter, 
and the magnetic moment of the neutron star. 

In §2 we set up the steady-state equations for a thin accretion disk around a fast X-ray pulsar and find analytic 
solutions. Torques on the neutron star due to accretion and magnetic drag on the disk are evaluated in §3. We 
summarize our results and draw conclusions in §4. 



2. MAGNETICALLY TORQUED ACCRETION 



2.1. Magnetospheric Accretion in the Thin Disk Model 



We now examine whether there exist solutions for a thin disk surrounding a "fast" pulsar, where the inner edge of 
the disk penetrates to the corotation radius r c . We start with the equation for conservation of angular momentum: 

where is the local Keplerian frequency in the disk, H is the full disk thickness, T a is the viscous torque per unit 
volume, and Tg is the magnetic torque on the disk per unit volume. In words, this equation says that the angular 
momentum transported inward by matter flow is driven by the viscous and magnetic torques, where a positive sign for 
r here indicates a tendency to slow down the disk rotation. In the standard "a-viscosity" prescription, the quantity 
T a is given by: 

'-■ikiw • < 4 > 
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where P is the pressure in the disk midplane; both P and H are functions of r, the radial distance from the neutron 
star. Here, aP is the standard viscous stress introduced by Shakura & Sunyaev (1974; hereafter SS). 

Before entering into the details of the model used below for the magnetic interaction between a disk and a mag- 
netosphere, consider first a simpler model (see also ST) which illustrates some properties of the problem that are 
independent of the detailed model used for the magnetic torques. 

In a standard thin disk accreting onto a slowly rotating object, the angular momentum flux, J, and the mass flux, 
M, are connected directly by the relation J — MrfOi, where r; is the radius of the inner edge of the disk and fij 
the Keplerian rotation frequency at this radius. This is independent of the detailed conditions at the inner edge 
of the disk, and is a consequence of the fact that the rotation curve has a maximum, near r;, where the viscous 
contrbution to the angular momentum flux vanishes (Lynden-Bell & Pringle 1972, SS). For accretion onto a rapidly 
spinning magnetosphere, there is no such extremum, and the angular momentum flux becomes an additional essential 
parameter of the problem. Its value is determined by the accretion rate and the details of the interaction between disk 
and magnetosphere, as illustrated in the following simplified problem. 



2.2. Accretion Torque Applied at Disk Inner Edge 



Consider the case of a rapidly rotating, magnetized neutron star that exerts a torque only over an infinitesimal radial 
range at the inner edge of the accretion disk - here taken to be located at r c . If we approximate the magnetic torque 
per unit volume as: 

T c S(r~r c ) 

= 2nrH (5) 
where 6(r — r c ) is a Dirac delta function and r c is the magnitude of the torque applied by the neutron star at the 
inner edge of the disk, we find that eq. (3) becomes: 

— - (Mtlr 2 ) = —- (PHr 2 ) + ^ ~ Tc) (6) 
2k Hr dr V ) Hr dr [ > + 2nrH ' W 

where a negative value of T c represents the torque from a "fast" pulsar. This equation can be directly integrated to 
yield: 

M (fir 2 - fl c r 2 c ) = 2-rraPHr 2 + T c , (7) 

where the evaluation of the PHr 2 term at r c is set equal to zero under the assumption that the density vanishes inside 
this point and therefore there can be no viscous torque acting on the innermost edge of the disk. We identify the 
combined terms (r c + MQ c r 2 ) as the net angular momentum flow into the neutron star, and define their sum to be 
jMtt c r 2 , where j is a dimcnsionless factor. The solution for the quantity PH can then be written as: 

PB-^U-lM ■ (8) 




Since we are considering only the case of fast pulsars, the factor j can, in principle, take on any value ^ 1 (i.e., T c ^ 0), 
the upper limit on j corresponding to zero magnetic torque on the inner edge of the disk. The physical meaning of 
this equation can be made clear without solving for the radial profiles of the individual parameters such as density, 
pressure, disk thickness, and so forth. This is facilitated by switching explicitly to a viscosity parameter rather than 
using the standard a prescription (SS). The relation between these two is: aP — —vprdQ/dr, where v is the coefficient 
of kinematic viscosity. With this substitution, and the explicit assumption of Keplerian rotation, eq. (8) becomes: 

where £ is the surface column density in the disk. Note, that in the transformation between v and a, they cannot 
both be taken to be constants independent of radial location, so that the radial profile of S would be different in the 
two formulations; nonetheless, the qualitative features that we wish to demonstrate will still be present in both. 

Figure n shows E(r), as a function of the parameter j. Note how the surface density near r c increases dramatically 
when the angular momentum flux deviates from its standard value of j = +1, and then drops abruptly to zero for 
r < r c . The inner regions of a disk with a net outward angular momentum flux (corresponding to j < 0) contain much 
more mass than in the standard case. For the net angular momentum flux (viscous plus advection) to be outward, the 
viscous flux has to be sufficiently large, which requires a high surface density (for further discussion see ST). 

These lowest-order considerations do not tell us, however, what determines the value of the angular momentum flux 
- except in the slow rotator case. A model for the transition region between disk and magnetosphere is needed, which 
specifies the magnetic torques on the disk in a physically plausible way. Such a model is developed below. One can 
already expect, however, that the equivalent of Fig. ^will look similar, but with the jump at r/r c = 1 broadened, 
while the parameter j will become a calculable function of the accretion rate. 
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Fig. 1. — Normalized surface density distribution S(r) in a standard thin disk with external torques acting only at the inner edge, 
i.e., at t = r c . Parameter j is the net specific angular momentum flux through the disk (advection plus magnetic torque) in units of 
the angular momentum at the inner edge. For accretion onto a slowly rotating object, j = +1 applies independent of conditions at the 
magnetoshere. For accretion onto a rapidly rotating magnetosphere, the angular momentum flux can be outward (j < 0, spindown of the 
star), and depends on details of the interaction between magnetosphere and disk. 



2.3. Model for a Disk with Distributed Magnetic Torques 

Here we take the magnetic torque on the disk to be distributed over a range of radial distances. Spcfically, we take 
the magnetic torque per unit volume to be: 

t b = . (10) 
2nH K ' 

Without a full, self-consistent magnetohydrodynamic solution of the disk equations, there is no simple way to determine 
B^. However, to a level of sophistication consistent with the SS-disk solutions, we simply adopt the following ad hoc 
but plausible prescription for B$: 

which we take to hold for r > r c , and oj s is the rotation frequency of the neutron star (see also, Wang 1987; 1995; 
1996 for more sophisticated and physically motivated versions of this relation). Here and in the following we assume, 
without loss of generality, that B z > 0. This expression has the property that the magnetic torque vanishes at the 
corotation radius, while B^ is comparable to B z at large distance. The latter is compatible with the fact that B^ can 
not be larger than B z over significant distances, for reasons of equilibrium and stability of the field above the disk 
plane. As the field is wound up by differential rotation between the star and the disk, the azimuthal component first 
increases, but the increasing energy in the azimuthal field component pushes the field configuration outward into an 
open configuration when the azimuthal component becomes comparable to the poloidal component. This was proven 
by Aly (1984, 1985) in a rather general context, and worked out in some detail for the case of a disk around a magnetic 
star by Lynden-Bell & Boily(1994). Finally, equation (3) reduces to: 

-^*(Mn 1 A="*{pH^)-^(l-Sl/v.) . (12) 
2-KHrdr\ J Hr dr v ; 2ttH y 11 K ' 

If we now integrate both sides of equation (12) with respect to r, starting at the corotation radius, we find: 

M (fir 2 - il c r 2 c ) = 2naPHr 2 - f £ (l - — ) dr , (13) 
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where we have taken the quantity PH to vanish at r c under the assumption that the magnetic field quickly sweeps 
up any matter deposited inside r c , thus maintaining a low pressure P. For the z component of the magnetic field we 
take p/r 3 to represent a dipole field whose axis is aligned approximately along the spin axis of the neutron star and 
perpendicular to the accretion disk. We further assume that f2 for all r > r c is always equal to the local Keplerian 
angular frequency. 

Before attempting to use equation (13) to solve for the disk structure, let us reflect a bit on its physical interpretation. 
We start by rewriting equation (13) schematically as follows: 

4W-4(» , o) = r«(r)-r B (r) , (14) 

where the four terms are as follows: J^i 7 ") 1S the angular momentum flux carried by matter across a boundary at 
radius r; Jj^{r c ) is the angular momentum flux carried by matter through the inner edge of the disk and onto the 
neutron star; —T a {r) is the net viscous torque on the disk interior to radius r; and T B (r) is the net magnetic torque 
on the disk interior to r. If we move the Ts(r) term to the other side of the equation and identify it as the angular 
momentum flux into the disk from the magnetic field of the neutron star, Jb, we find: 

T a {r) = j lk {r)-j lk {r c ) + j B {r) . (15) 

Written in this form, we can see that the viscous torques must transport angular momentum to the outer parts of the 
accretion disk at a rate equal to that transported by matter through r, minus the rate transported by matter through 
r c , plus the rate flowing out of the neutron star and into the disk (inside r) via the magnetic field. As r approaches 
r c the magnetic torque vanishes and J^{r) — > Jm{t c ), thereby requiring zero viscous torques at the inner boundary 
(P — > as well). At the outer edge of the disk, presumably tidal torques from the companion star remove the angular 
momentum flow driven by r Q , and feed it back into the orbit. The net angular momentum flow into the neutron star 
is Jj^(r c ) — Jb{t)- When this quantity becomes negative, the neutron star will be spun down. The larger the magnetic 
field, the larger must be the viscous torques, 2iraPHr 2 , in order to extract the increased angular momentum being 
injected into the disk. This is accomplished by the disk readjusting its structure, i.e., the radial profiles of P and H 
must change in such a way as to transport the additional angular momentum toward the outer parts of the disk. We 
show below that the principal enhancer of the viscous torque is a substantial increase in the pressure, i.e., the material 
becomes "dammed up" in the disk behind the corotation radius. 

If we now carry out the integral representing the magnetic torque on the disk (eq. 13), and solve for the quantity 
PH, wc find: 

PH = ^xF(r;r c ,p,M) , (16) 

where F is a dimensionless function of its argument and incorporates the effects of the magnetic torque on the disk. 
Specifically, F is: 



F(r;r c ,p,M) = 




where the first term in parentheses yields the SS solution in the absence of a magnetic field. The parameter £ is just 
the ratio of r m /r c . 

Three of the four remaining SS-disk equations remain unchanged. For the equation of state we take: P = pkT/m, 
where m is the mean molecular weight per particle. The vertical force balance equation is given by: P ~ GM pH 2 r~ 3 . 
Radiative transport in the vertical direction is represented by: T 4 ~ npHT*, where T and T e are the midplane and 
effective temperatures of the disk, respectively, and k is the radiative opacity evaluated at the disk midplane. Finally, 
we retain the fifth SS-disk equation for the heat dissipation per unit surface area of the disk: Q ~ aHPVl = uT^. 
In doing so, we have explicitly neglected the heat dissipation associated with the magnetic torque. Since, in the disk 
solutions we are proposing, the viscous torques must overcome the magnetic torques (which tend to repel the disk 
material), the concomitant viscous heating term should dominate that due to magnetic heating. For this reason, and in 
order to keep the disk equations in strictly algebraically solvable form, we neglect the magnetic heating. In all of these 
equations (except for the equation of state) we have neglected dimensionless coefficients of order unity. These factors 
have no impact on the form of the solutions for P, p, H, and T, and only a minor effect on the leading coefficients to 
the solutions. 

The four equations listed above, plus the equation for PH given in equation (16) can be solved algebraically in a 
manner analogous to that done for the original SS-disk equations to yield P, p, H, and T as functions of r, M, a, and 
p. In solving these equations we have taken k to be given by Kramers opacity (k~ 6 x 10 22 pT~ 3 ' 5 cm 2 gm _1 ) which 
is appropriate for most of the physical conditions found in our disk models, and the mass of the neutron star is fixed 
at 1.4 M Q . The solutions arc: 

P ~ 2 x 10 5 a - 9 / w Ml 7 6 /20 r- 21/8 F 17 / 20 dynes cm" 2 (18) 
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Fig. 2. — The density p as a function of radial distance, expressed in units of r c . Curves are for different values of the assumed accretion 
rate M, as indicated by the set of labels; the highest curve is for the highest value of M . 



/f~lx W 8 ^ 1 / 10 M^°r% S F 3 / 20 cm (19) 
T ~ 2 x lO^-^A^Vo 374 ^ 3710 K (20) 

p ~ 7 x lO-^" 7 / 10 ^ 1720 ^ 1578 ^ 11720 g cm" 3 (21) 

where Mie is the mass accretion rate in units of 10 16 gm sec -1 , and no is the radial distance in units of 10 10 cm. 
Note that for a vanishing magnetic held, the function F reduces to the function / 4 = 1 — y/r c /r, found in the SS-disk 
solutions (see eq. [17]). 

As an example of how the magnetic torque on the disk affects the disk solutions, we show in Figure [21 plots of 
p(r) for several different values of M. The plots are for the illustrative parameter values of a surface magnetic field 

B = 1 x 10 9 G, viscosity parameter a = 0.01, and a spin period P = 0.003 sec. The highest curve is for the value of M 
that corresponds approximately to the transition from a "slow" to "fast" pulsar; the lower curves are for progressively 
lower values of M and correspondingly "faster pulsars". Note that for r/r c ^> 1 the curves approach the standard 
SS solution with p cx M n /20 r -i6/s ( S ee eq.[21]). However, nearer the corotation radius, the density for the lower M 
values rises more steeply and is systematically larger than it would be in the absence of the magnetic field. 

Since the one disk equation that we have modified (from the original SS-formulation) involves the term PHr 2 (see 
eq. [12]), we can gain a better understanding of the solutions by plotting this quantity (see FigureOJ. Note that the 
behavior of the PHr 2 curves shows a flattening for lower values of M. For low values of M what happens is that the 
viscous torques on the disk (1st term on the r.h.s. of eq. [13]) become relatively large and go almost exclusively to 
counteracting the magnetic torques (tending to drive the disk material outward; 2nd term on the r.h.s. of eq. [13]). 
The difference between these two terms is quite small and provides a net drag on the disk material that drives the 
mass flow inward at a rate of M (l.h.s. of eq. [13]). In this limit (of low M), the £ 7 / 2 term in eq. (17) dominates 
over the (1 — J r c /r) term, and PHr 2 oc {1 — 3(r c /r) 3 + 2(r c /r) 9 / 2 }. For r only modestly larger than r c , PHr 2 — > 
constant, instead of growing approximately as yfr as in the slower pulsar case. These effects are clearly seen in the 
plot of PHr 2 vs. r shown in Fig. [3| 
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Fig. 3. — The product PHr 2 as a function of radial distance, 
expressed in units of r c . Curves are for different values of the 
assumed accretion rate M, as indicated by the labels. The curve 
for the highest (lowest) value of M is at the top (bottom). 
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Fig. 4. — The quantity 2naPH/UM as a function of radial 
distance, expressed in units of r c . This dimensionless ratio would 
equal f 4 = 1 — \Jr c /r and be independent of M for the simple case 
of an SS disk. Note that in this figure, the curve for the lowest 
(highest) value of M is at the top (bottom), i.e., the reverse of the 
order in Fig. [3] Values of the relative angular momentum flux into 
the neutron star j = J /(MQr 2 ) c are also indicated on the lower 
right. Negative values correspond to outward flux (spindown of 
the star) and high surface densities at the inner disk. 



As a further comparison of these magnetically torqued disk solutions against the SS solutions, we plot in Fig. |3]the 
quantity 2naPH/(flM) (~ 3-kvYi/M) which is a simple function of r for an SS disk (/ 4 = 1 — y/r c /r), approaches 
unity as r > r c , and is independent of M. Note that the curves which deviate most from unity (for intermediate and 
large r) are the solutions corresponding to the lower values of M, i.e., the faster pulsars. 

3. TORQUES ON THE NEUTRON STAR 

If, for fast X-ray pulsars, there exist accretion disks terminating at the corotation radius, as described above, then 
the net torque on the neutron star is straightforward to compute: 

r NS ~ My/ GMr c - jf B 2 z (r)r 2 (l - dr (22) 

(see also, LPP, GL, Wang 1987; 1995; 1996). Formally, the integral should be carried out from r c to an upper limit 
equal to the radius of the speed-of-light cylinder, r cy i. However, for the msec pulsars of interest r cy i S> r c . Equation 
(22) can be integrated to yield the torque on a fast X-ray pulsar of: 

r NS ~ My/ GMr c - £s . (23) 

Note that the first term on the right represents the usual spin-up torque due to matter accreting from r c , while the 
second term is a spin-down torque due to the magnetic field drag on the accretion disk. The factor of 1/9 in the 2nd 
term on the r.h.s. of eq. (23) results from our particular choice for (see, eq. [11]); for other plausible functional 
forms this factor would change somewhat. With our choice for B^, the torque reverses sign when r m — 9 2 ^ 7 r c ~ 1.87r c ; 
thus spinup can occur even while the pulsar is still 'fast', i.e., r rn > r c . In this regard, we remind the reader that 
the expression for r m given by eq. (1) is a formal definition of the magnetospheric radius and is not calculated self- 
consistently. 

In order to examine the behavior of the overall torque, Tns, as a function of M, we need a prescription for the 
case where the pulsar switches from "fast" to "slow", i.e., where r rn becomes smaller than r c . When this occurs we 
simply adopt the value of r m given by eq. (1) as the approximate location of the inner edge of the disk, and integrate 
the magnetic torque from r m to oo. By analogy with the "ad hoc" choice for B^ that we made for r > r c , we take 
i?0 = B z (l — u) s /Q) for r m < r < r c . Thus, for a slow pulsar, we find the torque on the neutron star to be: 



r NS ^ My/GMr, 1 



3r 



2_ 2 /r^\ 3/2 (rc 

3 V T rn J V T rTi 



(24) 




In Figure l^we show the results of evaluating \P/P\ as a function of M for a range of surface magnetic fields. In 
order to convert the torque, given in eqs. (23) and (24), to P/P one needs the rotation frequency of the neutron star. 
For this we have chosen an illustrative spin period of 2.5 msec, i.e., the value for the SAX 1808-3658 pulsar. The 
cusp-like structure in each curve is the location where P changes sign. The vertical bars show the approximate range 
of M inferred for SAX 1808-3658 during its most recent outburst in 2002 (for a distance of ~ 4 kpc; Chakrabarty 
& Morgan 1998; Homer et al. 2001). The heavy arrow is the average value of P/P ~ 2 x 10~ 8 yr _1 observed over 
a period of ~4 weeks during the 2002-outburst while the source declined in intensity by a factor of ~50 (Morgan, 
Galloway, & Chakrabarty 2003). By inspection of Figure [S] it can be seen that the proposed model indicates a surface 
B field of ~ 8 x 10 8 G. We point out that the magnitude of the observed torque is substantially larger than that 
expected for the matter torque alone over most of the duration of the outburst cycle (especially for the lower values 

of M). In this respect, the proposed model can successfully explain the observed spin period behavior - especially if 
the magnetic field is ultimately confirmed to be as large as ~ 8 x 10 8 G. We also note, however, that a reanalysis of 
the spin behavior of SAX 1808-3658 during its 1998 outburst (Morgan, Galloway, & Chakrabarty 2003) has detected a 
significant spinup over most of its on state — and of a magnitude that is even larger than that of the recently observed 
spindown. Thus, we caution that the msec pulsars are apparently capable of exhibiting spinups and spindowns that 
are inconsistent with the inferred value of M, as is the case for their more slowly rotating cousins (see, e.g., Nelson et 
al. 1997). 

XTE J0929-314 has thus far been observed only during a single outburst (Galloway et al. 2002). During that time 
the source underwent spindown with a value of P/P ~ 1.6 x 10 -8 yr _1 , and had a luminosity L x < 1.5 x 10 36 ergs s _1 
(for an assumed distance of 5 kpc; Galloway et al. 2002). An analysis similar to that done for SAX 1808-3658 (Fig. 
|SJ) yields a comparable value for the magnetic field of ~ 1.0 x 10 9 G. Accretion torques have not yet been reported for 
the other three known msec X-ray pulsars. 
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4. DISCUSSION AND CONCLUSIONS 

In this work, we have explored the proposition that fast X-ray pulsars can continue to accrete from a thin disk, even 
for accretion rates that place the nominal magnetospheric radius well beyond the corotation radius. We hypothesize 
that the inner edge of the accretion disk will be located just inside the corotation radius. We modify the Shakura- 
Sunyaev disk equations by adding a simple magnetic torque prescription, and find analytic solutions to these equations. 
The density and pressure profiles are significantly modified over the SS solutions as the fastness parameter increases. 
We have not analyzed the physical stability of these solutions, and therefore cannot comment on the lower limit 
on M which would be allowed in our model. The form of the magnetic torques we have assumed are admittedly 
subject to substantial uncertainty; however, the main conclusions we reach, e.g., accretion over a wide range of M and 
compatibility of spindown during accretion should be largely independent of the exact form adopted (as argued before 
in ST). 

The proposed model directly confronts the observational fact that the accretion by millisecond X-ray pulsars ap- 
parently persists over a wide range of M, including down to very low values of M < lO~ n M0 yr" 1 . The model can 
also account for the magnitude of the observed spindown torques that have been observed for two of the msec X-ray 
pulsars. It does not easily explain the magnitude of the spinup episode in SAX J1808-3658 during its 1998 outburst. 
The observed spindown behavior in SAX J1808-3658 (during the 2002 outburst) and XTE J0929-314 requires surface 
magnetic fields of ~ 8 x 10 8 G; this is a prediction that can possibly be used to verify or falsify the model. Also, if 
correct, our model would remove the constraints on the neutron star radius in SAX J1808-3658 found by Burderi & 
King (1998) and Psaltis & Chakrabarty (1999). 

Campana et al. (2001) reported a rapid transition in the X-ray luminosity of 4U 0115+63 by a factor of > 250 over 
a 15-hr interval, during which time the pulsations (P = 3.6 sec) continued with only minor changes in pulse fraction. 
This change in luminosity, according to the standard accretion picture, should have signaled the onset of a transition 
from direct accretion onto the neutron star surface to a propeller mode. Campana et al. (2001) developed a simple 
model to account for this behavior without invoking the type of "fast-pulsar accretion" advocated in this work. In 
their model, the actual value of M approaching r m changes by only a factor of a few, while the luminosity changes by 
a factor of > 250. Based on the results presented, it is difficult to judge whether this model is to be preferred over the 
one presented herein. First, we note that these authors did not allow for the possibility of direct, continuous accretion 
onto a fast pulsar. Second, they assert that the rapid change in luminosity over a 15-hr interval cannot naturally be 
accounted for by variations in M flowing through the disk. However, given our lack of detailed understanding of Be 
star excretion disks, it does not seem completely implausible that the change in X-ray luminosity directly reflects the 
change in M into the disk. Finally, interesting predictions of both models would be the sign and magnitude of the 
spin-torque behavior of the neutron star — however, no such measurements were reported. 

Twin-peak kHz quasiperiodic oscillations ("QPOs") have been observed in the X-ray intensity of two of the msec 
pulsars (SAX 1808-3658 and XTE J1807-294; Chakrabarty et al. 2003; Markwardt et al. 2003). In some models, 
one of the two kHz QPO peaks represents the orbital frequency of blobs of matter near the neutron star (see, e.g., 
Miller, Lamb, & Psaltis 1998; Stella & Vietri 1999). We note that such high Keplerian frequencies cannot be reached 
outside of r c where the accretion disk is located in our model. If such an association between the kHz QPOs and 
Keplerian motion is firmly established, our model would not naturally accommodate these orbital frequencies; however, 
we emphasize that the kHz QPOs are not well understood at this time. 

At an accretion rate of M = 2 x 10 15 gm s" 1 ~ 3 x 10~ 11 M Q yr" 1 our disk model requires <~ 2 yr to fill from an 
empty state; the corresponding filling time of an SS disk is <~ 1 yr. The difference of <~ 1 yr is comparable to the 
recurrence time of these transient msec pulsars, and may in some way be related to the transient outbursts. This 
aspect of the model should be pursued with time-dependent calculations. 



We thank Alan Levine and Eric Pfahl for helpful discussions. Edward Morgan and Deepto Chakrabarty supplied 
important results on the spin behavior of SAX 1808-3658 prior to their publication. James Kiger and John Belcher 
participated in an earlier version of this research. One of us (SR) acknowledges support from NASA ATP Grant 
NAG5-12522. 
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